Problem Note 60270: PROC GLIMMIX might compute incorrect least squares means when an EFFECT statement is used and the AT= option is specified in the LSMEANS statement
 |  |  |  |
PROC GLIMMIX produces incorrect least squares means when
- the MODEL statement includes a constructed effect defined in an EFFECT statement,
- more than one effect is specified in an LSMEANS statement, and
- the AT option is specified in the LSMEANS statement
The following code illustrates the combination of statements and options that cause the problem:
proc glimmix;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
lsmeans A|B / at x=10;
run;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
lsmeans A|B / at x=10;
run;
To circumvent the problem, specify a single effect in the LSMEANS statement when using the AT= option:
proc glimmix;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
lsmeans A / at x=10;
lsmeans B / at x=10;
lsmeans A*B / at x=10;
run;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
lsmeans A / at x=10;
lsmeans B / at x=10;
lsmeans A*B / at x=10;
run;
Another circumvention is to use the PLM procedure for computing the LSMEANS:
proc glimmix;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
store sasuser.glimmix;
run;
proc plm restore=sasuser.glimmix;
lsmeans A|B / at x=10;
run;
class A B sub;
effect spl = spline(x);
model y = A B A*spl;
random A*B / subject=sub;
store sasuser.glimmix;
run;
proc plm restore=sasuser.glimmix;
lsmeans A|B / at x=10;
run;
Operating System and Release Information
| Product Family | Product | System | Product Release | SAS Release | ||
| Reported | Fixed* | Reported | Fixed* | |||
| SAS System | SAS/STAT | z/OS | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 |
| z/OS 64-bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft® Windows® for x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8 Enterprise 32-bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8 Enterprise x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8 Pro 32-bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8 Pro x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8.1 Enterprise 32-bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8.1 Enterprise x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8.1 Pro 32-bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 8.1 Pro x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows 10 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows Server 2008 | 14.1 | 9.4 TS1M3 | ||||
| Microsoft Windows Server 2008 R2 | 14.1 | 9.4 TS1M3 | ||||
| Microsoft Windows Server 2008 for x64 | 14.1 | 9.4 TS1M3 | ||||
| Microsoft Windows Server 2012 Datacenter | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows Server 2012 R2 Datacenter | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows Server 2012 R2 Std | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Microsoft Windows Server 2012 Std | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Enterprise 32 bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Enterprise x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Home Premium 32 bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Home Premium x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Professional 32 bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Professional x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Ultimate 32 bit | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Windows 7 Ultimate x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| 64-bit Enabled AIX | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| 64-bit Enabled Solaris | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| HP-UX IPF | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Linux for x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
| Solaris for x64 | 14.1 | 14.3 | 9.4 TS1M3 | 9.4 TS1M5 | ||
Incorrect least squares means are generated in the GLIMMIX procedure when the model includes a constructed effect and multiple LS-means are requested at a specific level of a covariate.
| Type: | Problem Note |
| Priority: | high |
| Topic: | Analytics ==> Mixed Models SAS Reference ==> Procedures ==> GLIMMIX |
| Date Modified: | 2017-08-28 11:45:50 |
| Date Created: | 2017-04-07 14:34:15 |